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METHOD  OF  ANALYSIS  FOR  COMPRESSIBLE  FLOW  PAST  ARBITRARY 
TURBOMACHINE  BLADES  ON  GENERAL  SURFACE  OF  REVOLUTION 
By  Chung-Hua  Wu  and  Curtis  A,  Brown 


SUMMARY 

A  method  of  obtaining  the  flow  of  a  nonviscous  coii5)ressible  fluid 
past  arbitrary  compressor  or  turbine  blades  between  two  neighboring  sur¬ 
faces  of  revolution  is  presented.  The  equations  of  continuity  and 
motion  obtained  for  such  flow  are  combined  into  a  nonlinear  second- 
order  differential  equation  in  terms  of  a  stream  function  defined  for 
such  flow.  Numerical  solution  by  the  use  of  differentiation  coeffi¬ 
cients  for  unequally  spaced  grid  points  is  suggested.  Means  for  satis¬ 
fying  the  boundary  conditions  outside  the  channel  and  solutions  by  the 
relaxation  method  with  manual  computation  and  by  a  matrix  method  with  a 
large-scale  digital  computing  machine  are  described. 

Satisfactory  results  were  obtained  by  the  method  described  in  the 
investigation  of  the  detailed  flow  variation  of  a  compressible  fluid 
past  typical  high-solidity^  highly  cambered  thick  turbine  blades  on  a 
cylindrical  surface.  The  variations  in  fluid  properties  across  the 
channel  appeared  to  be  representable  by  a  second-degree  function.  The 
mean  streamline  approximately  followed  the  shape  of  the  mean  channel 
line  of  the  cascade  and  had  lower  curvature.  The  variation  of  specific 
mass  flow  along  the  mean  streamline  followed  the  trend  of  the  varia¬ 
tion  in  the  channel  width.  In  general,  the  variation  in  specific  mass 
flow  was  significantly  higher  than  that  given  by  the  ratio  in  channel 
width  and  the  effect  extended  outside  the  channel.  The  velocity  distri¬ 
bution  around  the  blade  obtained  in  the  theoretical  calculation  compares 
very  well  with  experimental  values. 


INTRODUCTION 

A  basic  aerodynamic  problem  of  turbojet  and  turbine-propeller 
engines  is  the  flow  of  compressible  fluid  past  a  series  of  blades  in 
circular  arrangement.  In  axial-flow  type  turbomachines,  if  the  blades 
are  relatively  short  in  the  radial  dimension  and  are  bounded  by  cylin¬ 
drical  walls,  the  theoretical  flow  passing  through  the  blades  is  usually 
computed  on  the  basis  of  two-dimensional  flow  on  a  cylindrical  surface. 
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which  is  then  developed  into  a  plane  for  convenience  of  calculation.  A 
number  of  methods  have  been  proposed  to  obtain  the  theoretical  flow 
throu^  such  a  given  cascade  of  airfoils,  including  analytical  methods 
using  conformal  mapping,  interference  technique,  and  Prandtl-Glauert  or 
Chaplygin-K^mAn-Tsien  approximation  for  compressible  flow  (references  1 
to  6),  graphical  procedure  (references  7  to  lO),  and  other  mechanical 
and  electrical  devices  (reference  10).  For  radial-  and  mixed-flow  type 
turbomachines,  the  available  methods  consider  incompressible  flow  in  a 
radial  plane  (references  7  and  11  to  16)  or  compressible  flow  on  a  coni¬ 
cal  surface  (reference  17).  However,  in  current  axial -flow  turbo¬ 
machines,  the  hub  and  casing  walls  may  be  either  tapered  or  cinrved, 
which  causes  the  fluid  to  flow  on  a  noncylindrical  surface;  and  in  cen¬ 
trifugal  machines,  the  flow  surfaces  are  usually  quite  curved  toward 
the  inlet  (see  fig.  l) .  A  method  of  analysis  is  therefore  developed  at 
the  MCA  Lewis  laboratory  to  analyze  the  two-dimensional  compressible 
flow  for  the  fluid  between  two  neighboring  surfaces  of  revolution  in 
these  turbomachines  (fig.  l) .  The  equations  of  continuity  and  motion 
for  irrotational  absolute  flow  are  first  obtained  for  such  flows  and  are 
then  combined  into  a  nonlinear  second-order  partial  differential  equa¬ 
tion.  Because  the  change  of  fluid  properties  passing  through  turbo¬ 
machine  blades  is,  in  general,  large  and  the  shapes  of  surface  and 
blades  are  arbitrary,  numerical  solution  by  the  finite-difference 
approach  is  suggested. 

Solving  the  flow  through  a  cascade  of  blades  of  arbitrary  camber 
and  thickness  by  the  finite-difference  method  involves  two  main  diffi¬ 
culties:  (l)  the  curved  boundaries  formed  by  the  blade  surface,  and 

(2)  the  large  number  of  grid  points  necessary  to  cover  the  whole  flow 
region.  The  first  difficulty  is  removed  by  the  recently  available  dif¬ 
ferentiation  formula  and  coefficients  for  unequally  spaced  grid  points 
(reference  is)  •  With  these  fonnulas  and  coefficients,  the  curved  bound¬ 
ary  can  be  handled  in  the  same  manner  as  a  straight  bOTUidary.  The 
second  difficulty  can  be  reduced  by  the  use  of  higher-degree  polynomial 
representations  (references  18  and  19) .  Furthermore,  if  a  modern  high¬ 
speed  large-scale  digital  machine  is  available,  the  set  of  simultaneous 
difference  equations  in  the  stream  function  can  be  very  quickly  solved 
by  the  matrix  process  given  in  reference  18,  and  consequently,  suc¬ 
cessively  iniproved  solutions  for  compressible  flow  can  be  obtained  in  a 
reasonably  short  time.  Without  such  a  machine,  it  is  desirable  to 
obtain  an  approximate  solution  for  con5)ressible  flow  by  either  graphic, 
mechanical,  or  other  approximate  methods  and  to  use  the  relaxation 
method  (references  20  and  21)  for  the  final  Improvements. 

In  turbomachines  with  a  relatively  large  radial  dimension,  curved 
hub  and  casing  walls,  or  if  designed  on  velocity  diagrams  other  than  the 
free-vortex  type,  the  radial  flow  assumes  primary  importance.  For  such 
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turlDoma chines,  methods  of  flow  analysis  have  been  proposed  that  take 
into  account  the  radial  flow  but  consider  only  the  variation  of  the 
fluid  properties  on  a  mean  stream  surface  passing  between  two  blades 
(reference  22).  When  applying  such  methods  of  analysis  to  turbomachines 
with  thick  blades,  it  is  necessary  to  have  some  knowledge  of  the  effect 
of  area  reduction  due  to  blade  thickness  on  the  mean  flow  and  the  rela¬ 
tion  between  the  blade  shape  and  the  shape  of  the  mean  streamline  at 
various  radii.  For  this  reason,  the  method  developed  is  applied  to 
investigate  the  detailed  compressible  flow  in  a  typical  turbine  cascade 
with  highly  cambered  thick  blades.  The  theoretical  velocity  on  the 
blade  is  compared  with  experimental  data  obtained  at  the  NACA  Lewis 
laboratory. 


SYMBOLS 

The  following  symbols  are  used  in  this  report: 
a,b,c  points  on  streamline 
B  differentiation  coefficients 

C  constant 

E  error 

H  total  enthalpy  based  on  absolute  velocity 

h  static  enthalpy 

I  h  +  I  W2  -  i  (X)2r2 

L  blade  length  projected  on  z-axis 

7,cp  orthogonal  coordinates  on  mean  surface  of  revolution, (fig.  l) 

M  total  mass  flow  passing  throu^  space  between  two  neighboring 

blades 

n  degree  of  polynomial 

P  pitch  or  spacing 

p  static  pressure 

r  radial  distance  from  the  axis  of  machine^  (fig-  l) 
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t  “blade  thickness  in  pitch  direction 

W  velocity  relative  to  hlade 

y 

z  distance  along  axis  of  machine^  (fig*  l) 

1  1  Wcp 

a  flow  angle  in  Zcp-plane,  tan"-^-  -  ^ 

r  w-^ 

X  ratio  of  specific  heats 

6  grid  spacing 

p  mass  density 

1  r 

a  tan-J-  (fig*  1) 

T  normal  thickness  of  stream  filament  of  revolution,  (fig.  l) 

i|r  stream  function 

CD  angular  velocity  of  hlade 

Subscripts: 


i 


m 

n 

P 

s 

T 


inlet 

grid  points  where  differentiation  coefficients  are  applied 

meridional  and  circumferential  component,  (fig.  l) 

mean  streamline 

degree  of  polynomial 

pressure  surface 

suction  surface 

total,  or  stagnation,  state 
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Superscripts: 

a,lD,c,  .  .  .  grid  points 

1  2  indicate  first  and  seccsnd  derivatives,  respectively 


METHOD  OF  SOLUTION 

The  steady  isentropic  flow  of  a  nonviscous  compressible  fluid  along 
a  stream  filament  of  revolution  described  by  the  orthogonal  coordinates 
I  and  cp  on  the  mean  surface  of  revolution  and  a  varying  normal  thick¬ 
ness  of  the  filament  T  (fig.  l)  is  governed  by  the  following  equa¬ 
tions  of  continuity  and  irrotational  absolute  flow  and  the  isentropic 
pressure-density  relation  (see  appendix  for  derivation  of  the  first  two 
equations) : 


a(TPW^r)  a(TpWcp)  ^ 

—  +  '  sip  ^  ° 


SWq,  1  SWj  VUsin  cT 

Sr  "  r  - +  3co  sin  a  =  0  (2) 
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In  the  nmerical  solution  of  the  problem,  iteration  is  necessary 
for  compressible  flow.  A  network  of  grid  points  is  chosen  to  cover  the 
whole  flow  region,  and  finite-difference  equations  are  obtained  for  the 
first  three  terms  of  equation  (s) .  The  last  term  is  taken  as  constant 
during  the  solution  of  the  t  values  and  recomputed  from  the  improved 
values  of  \lr  for  the  next  solution  (references  20  and  2l) .  In  this 
way,  the  numerical  coefficients  at  each  grid  point  are  determined  by 
the  geometry  of  the  problem  alone. 

The  value  of  density  is  most  conveniently  expressed  throughout  the 
flow  region  in  terms  of  velocity  and  enthalpy  by  the  use  of  equation  (3) 
and  the  relation  for  a  perfect  gas: 


1 


where  I  is  constant  throughout  the  flow  region  for  adiabatic,  abso¬ 
lutely  Irrotational,  and  steady  relative  flow  (reference  22) .  (For 
application  to  stator  blades,  o)  becomes  zero,  W  becomes  the  abso¬ 
lute  velocity,  and  I  becomes  H. )  If  the  conponents  of  W  as 
expressed  by  equations  (4)  are  substituted  in  equation  (6),  the  follow¬ 
ing  equation  is 


-P 

PT,1 


I 


+  i 

“hI 


2  Hi 


_1 


(7) 


In  order  to  facilitate  the  evaluation  of  the  density  from  the 
t-derivatives,  equation  (?)  can  be  rewritten  in  the  following  form: 


2 

Z  =  (1  - 


(7a) 


where 


(7b) 


2151 


2151 


MCA  TN  2407 


7 


and 


0  = 


^-1 


2(t  PT,i)^  % 


y+1 

r~i 


With  the  relation  between  p  and  ^jr-derivatives  in  this  form^  either  a 
table  may  be  constructed  by  computing  S  for  equally  spaced  values  of  $ 
through  Iteration  on  a  punch  card  machine^  or  a  graph  may  be  constructed 
by  computing  $  for  a  nimber  of  S.  Because 

and  T  are  all  known  values  ^  with  the  aid  of  the  table  or  graphs  the 
value  of  density  at  any  point  can  be  easily  obtained  after  the  first- 
order  derivatives  of  ^  are  computed  by  using  the  differential  coeffi¬ 
cients  for  equally  and  unequally  spaced  grid  points  given  in  refer¬ 
ences  23  and  18^  respectively. 


The  first-order  derivatives  of  density  are  then  computed  and  com¬ 
bined  with  the  first-order  derivatives  of  If  term  in  equa¬ 

tion  (5)  is  denoted  by  J  and  the  differentiation  coefficients^  which 
multiply  the  \lr  value  at  point  j  to  give  the  second  and  first 
derivatives  at  point  i  using  an  nth-degree  polynomial^  are  denoted  by 

and  ^B^^  respectively^  the  finite-difference  form  of  equation 

n  j  ^  j 

at  the  grid  point  i  can  be  written  as 


n 


j=0 


Li  ^  [sin_a 
|n  j  \  r 


5  la 

~§r7 


J 


n 


k=0 


nft 

(^)2 


k 


+ 


Ji  =  0 


where  and  denote  ^  values  on  the  surface  along  I  and  cp 

coordinates,  respectively  (figs.  1  and  2).  If  second-degree  polynomial 
representation  is  chosen  (n=2),  three  points  are  involved  in  the 
numerical  differentiation  in  each  direction,  and  the  center-point  for¬ 
mula  can  always  be  used  for  all  interior  grid  points.  If  fourth-degree 
polynomial  representation  is  chosen  (n=4),  five  points  are  involved 
in  each  direction  and  the  off-center-point  formula  has  to  be  used  at 
the  point  next  to  the  boundary.  Two  typical  cases  where  odd  spacings 
are  used  in  the  calculation  of  the  coefficients  because  of  a  curved 
boundary  are  shown  in  figure  2.  In  the  case  where  ilfl  is  at  j=2, 
k=2,  the  center-point  formula  is  used  in  both  directions,  with  equal- 
space  7 -direction  coefficients  obtainable  from  reference  23  and 
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with  cp-directlon  coefTiclents  for  the  odd  space  at  the  boundary  obtain¬ 
able  from  reference  18.  In  the  other  case  where  ilr^  is  at  j=2, 
k=3,  again  the  equal-space  center-point  coefficients  are  used  in  the 
i -direction;  but  in  the  q? -direction  an  off -center-point  formula  for 
coefficients  must  be  used,  including  an  odd  space  at  the  boundary, 
again  obtainable  from  reference  18.  It  is  desirable,  where  possible,  to 
have  the  spacing  near  the  boundary  somewhat  shorter  than  the  regular 
equal  spacing  so  that  accuracy  at  the  point  next  to  the  boundary  is 
comparable  with  that  at  other  points.  By  using  a  five-point  system, 
many  less  grid  points  are  required  in  comparison  with  the  three-point 
system  for  the  same  accuracy.  (See  reference  18.)  In  the  present 
problem,  the  variation  of  ijr  is  only  rapid  near  the  leading  and  trail¬ 
ing  edges  of  the  blade  and  the  grid  pattern  should  be  determined  for 
these  regions  first. 

After  the  degree  of  polynomial  representation  and  the  size  of  grid 
have  been  chosen,  a  finite -difference  equation  (8)  can  be  easily  written 
for  each  interior  grid  point  using  the  differentiation  coefficients  for 
equal  and  .’neqioal  grid  spacings.  The  constant  t  value  on  the  suction 
surface  can  be  arbitrarily  chosen  and  the  constant  \lf  value  on  the 
pressure  surface  is  determined  by  the  mass  flow  passing  through  the 
channel  or  at  the  inlet  as  follows: 


- 


^s  = 


‘<Pp 


Si|r  , 


(ft 


TpW^r  dcp  =  M 


(9) 


Outside  the  channel  two  reference  lines,  a  pitch  angle  apart,  can  be 
drawn  either  parallel  to  the  Z-axis  or  parallel  to  the  inlet  and  exit 
angles  in  the  2cp-plane.  (Parallelism  with  the  inlet  and  exit  angles, 
as  shown  in  fig.  3,  is  a  little  more  convenient  for  the  purpose  of 
drawing  streamlines  after  the  solution  is  obtained. )  Because  the 
fluid  state  repeats  itself  for  each  pitch  angle  along  the  Cp-dlrection, 
only  the  grid  points  lying  between  the  reference  lines  need  to  be 
included  in  the  calculation  when  the  number  of  grid  spaces  per  pitch 
angle  is  an  integer,  and  the  central  five-point  formula  can  always  be 
used  at  every  grid  point  between  these  lines.  For  example,  at  grid 
point  c  (fig.  3),  the  vlr  values  at  a  and  b  are  obtained  in  terms 
of  the  ijf  values  at  1  and  j,  which  are  a  pitch  angle  apart  from 
a  and  b,  respectively,  by  the  following  relation: 


¥ 


^l^a  =  itfi  +  (ifp  -  ts) 

i|,b  = 


(10) 


ir 
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Inlet  and  Exit  Boundary  Conditions 

The  solution  of  the  present  problem  in  terms  of  t  has  an  inter¬ 
esting  feature  in  the  boundary  conditions.  As  noted  in  the  previous 
section,  the  value  is  constant  along  either  the  suction  or  pressure 
surfaces,  which  are  fixed  boundaries^  and  outside  the  channel,  there  is 
no  longer  a  fixed  boundary,  but  there  is  the  condition  that  the  flow 
repeats  itself  for  every  pitch  angle  in  the  cp -direction.  Sufficiently 
far  upstream  of  the  blade,  an  inlet  station  can  be  chosen  so  the  bound¬ 
ary  condition  there  requires  the  fluid  state  to  be  uniform  in  the 
cp-direction  and  prescribes  the  inlet  angle  Pj_.  The  following  two 
methods  have  been  devised  to  account  for  this  boundary  condition. 

Fixed-angle  method.  -  A  part  of  the  grid  system  in  the  Icp-plane 
near  the  inlet  station  1-i  is  shown  in  figure  4(a) .  In  order  to  write 
the  finite -difference  equation  at  point  a,  the  value  of  t  at 
point  c  to  the  left  of  station  i-i  can  be  obtained  by  using  the  given 
inlet  angle  pp  and  the  linear  variation  of  ilf  at  station  i-i 

67  1:311  CL  A 

^lfC  =  ^d  =  ^a  +  - £  (.[fb  _  ^a) 

Ocp 


The  coefficients  of  \lf®'  and  in  equation  (ll)  after  having  been 

multiplied  by  are  added  to  the  regular  coefficients  at  points  a 

and  b;  otherwise  the  points  a  and  b  are  treated  in  the  same  manner 
as  other  interior  points.  Whether  the  inlet  station  1-i  is  chosen  suf¬ 
ficiently  far  away  from  the  blade  is  indicated  by  the  linearity  of  the 
t  variation  at  that  station  as  obtained  in  the  solution. 

Streamline -adjustment  method.  -  If  the  first  station  i-i  is  chosen 
sufficiently  far  from  the  blades,  the  variation  of  stream  function  to 
the  left  of  the  station  i-i  is  linear  in  the  circumferential  direction. 
The  value  of  the  stream  function,  however,  depends  on  the  inlet  angle. 

If  solutions  for  a  range  of  inlet  angles  are  desired,  they  can  be 
obtained  by  specifying  a  number  of  linearly  varying  stream  functions  to 
the  left  of  station  i-i  as  fixed  boimdary  values.  The  streamline 
obtained  in  the  solution  then  gives  the  value  of  the  inlet  angle.  If, 
however,  the  solution  for  a  certain  inlet  angle  is  desired,  the  stream¬ 
line  obtained  in  the  solution  must  be  adjusted  according  to  that  inlet 
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angle,  (for  example  as  ab  in  fig.  4(1)  is  adjusted  to  position  ac), 
thereby  obtaining  an  improved  set  of  boundary  values  of  the  stream 
functions  to  be  used  in  the  next  calculation.  This  method  is  of  course 
not  as  accurate  and  convenient  as  the  fixed-angle  method  for  obtaining 
a  solution  for  a  given  inlet  angle  but  is  desirable  in  the  matrix  solu¬ 
tion  (to  be  discussed  in  the  following  section)  because  the  inlet  angle 
is  not  involved  in  the  matrix  factorization,  making  the  same  matrix 
factors  usable  for  a  range  of  inlet  angles  and  Mach  numbers. 

At  the  exit  station  far  downstream  of  the  blade,  the  same  methods 
can  be  applied.  For  a  blade  having  a  sharp  trailing  edge,  the  Kutta- 
Joukowski  condition  can  be  used  and  the  correct  exit  angle  far  down¬ 
stream  is  the  one  that  gives  the  flow  at  the  trailing  edge  satisfying 
that  condition.  For  round  trailing  edges,  either  the  position  of  the 
stagnation  point  can  be  assumed  or  some  available  empirical  rules  for 
the  exit  angle  used.  If  the  calculation  is  made  to  compare  with  cer¬ 
tain  experimental  results,  the  measured  exit  angle  can  be  used. 


Solution  of  Finite-Difference  Equations 

With  the  grid  system  and  the  degree  of  polynomial  representation 
chosen  and  the  boundary  conditions  taken  into  account,  the  problem 
remaining  is  the  solution  of  the  set  of  N  linear  algebraic  equa¬ 
tions  (8)  written  for  N  interior  grid  points.  For  a  small  number  of 
solutions  with  a  given  blade,  the  best  method  is  the  relaxation  method 
(reference  20).  A  modification  of  this  method  involving  the  use  of 
higher-order  differences  is  suggested  by  Fox  (reference  18).  Formulas 
and  the  table  of  coefficients  obtained  in  reference  18  enable  the 
direct  use  of  higher-degree  polynomials  for  problems  with  curved  bound¬ 
aries.  For  the  present  flow  problems,  it  is  necessary  to  include  a 
large  domain  to  get  to  the  boundary  conditions  that  are  given  at  places 
far  from  the  blades.  The  use  of  higher-degree  polynomials  will  greatly 
reduce  the  numerical  work. 

If  a  number  of  cases  are  to  be  solved  for  a  given  set  of  blades 
on  a  given  surface,  it  is  advantageous  to  solve  the  problem  on  a  large- 
scale  digital  conrputing  machine.  If  a  very  hi^- speed  digital  machine 
is  available,  the  simultaneous  equation  is  best  solved  by  Liebmann's 
iterative  process,  which  is  the  most  simple  to  set  up  (reference  24). 

If  only  a  relatively  slow-speed  machine  is  available,  the  matrix 
process  suggested  in  reference  18  is  most  suitable. 
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Matrix  solution.  -  In  the  matrix  solution,  the  entire  equation  (S) 
is  used  with  p  being  kept  constant  during  any  one  cycle.  Several 
cycles  are  necessary  for  compressible  flow  because  p  is  in  a  continual 
process  of  change  during  each  calculation  according  to  equation  (7). 

The  matrix  method  is  found  most  satisfactory  when  several  solutions  are 
required  for  a  given  blading  configuration,  (for  instance,  in  the  case 
where  the  flow  is  to  be  determined  throu^  a  blading  at  several  inlet 
Mach  numbers  or  inlet  angles)  because  the  factorization  of  the  matrix 
is  the  most  time-consuming  part  of  the  problem  and  once  completed  the 
matrix  may  be  used  again  and  again  as  one  of  the  inlet  flow  conditions 
is  varied.  For  example,  in  a  problem  similar  to  the  cascade  problem 
in  which  400  grid  points  and  foTirth-degree  polynomials  were  used,  it 
took  about  60  hours  to  factorize  the  matrix  on  an  IBM  Card  Programmed 
Electronic  Calculator,  and  solving  the  set  of  values  of  ilr  at  all 
grid  points  for  a  given  set  of  J  values  only  takes  2  hours.  The 
rounding-off  error  is  found  to  be  very  small. 

Relaxation  solution.  -  In  the  relaxation  solution,  the  entire  equa¬ 
tion  (s)  is  used  in  the  calculation  of  the  residuals}  but  only  the  first 
two  terms  are  used  in  the  actual  relaxation.  As  in  the  matrix  solution, 
several  cycles  are  required  for  compressible  flow  because  p  is  not 
allowed  to  vary  dirring  the  relaxation  process  whereas  actually  p  is 
changing  according  to  equation  (7).  Consequently,  after  each  cycle 
new  values  of  p  are  calculated  and  used  in  the  determination  of  an 
improved  J  as  in  the  matrix  method. 

Iteration  method  of  Liebmann.  -  The  iteration  method  of  Liebmann 
is  simplest  to  set  up  on  a  digital  computing  machine;  but  it  is  the 
slowest  method,  and  therefore  it  can  only  be  used  for  the  present 
problem  on  a  high-speed  machine.  The  boundary  conditions  and  the 
between-cycles  J  correction  process  are  the  same  as  those  in  the 
matrix  and  relaxation  methods. 


ACCURACY  OF  FIKAL  RESULTS 

Several  factors  affect  the  accuracy  of  the  final  results:  the 
degree  of  Lagrangian  polynomial  used,  the  size  relative  to  the  airfoil 
of  the  network  spacings,  the  smallness  to  which  the  residuals  are 
reduced,  and  the  number  of  cycles  carried  out  for  compressible  flow. 

The  fourth-degree  Lagrangian  polynomial  has  been  found  to  be  the 
most  applicable  to  the  type  of  calculation  under  consideration.  The 
use  of  a  second-degree  polynomial,  in  general,  requires  a  smaller 
grid  spacing,  necessitating  a  much  larger  number  of  network  points  than 
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the  fourth -degree  polynomial.  Polynomials  of  higher  degree  than  the 
fourth  are  not  desirable  within  the  channel,  because  fitting  the  points 
necessary  for  these  higher-degree  polynomials  circumferentially  Inside 
the  channel  may  require  a  smaller  grid  spacing.  Also  the  accuracy  near 
the  leading  and  trailing  edges  of  the  blade  Is  best  obtained  by  using 
RTna  1 1  spacing.  Outside  the  channel,  however,  the  use  of  a  higher- 
degree  polynomial  would  not  only  reduce  the  amount  of  work  but  also 
would  specify  more  conipletely  the  condition  that  the  flow  repeats  cir¬ 
cumferentially  every  pitch  angle. 

Formulas  that  can  be  used,  after  the  problem  has  been  completed, 
to  determine  the  order  of  magnitude  of  the  error  at  any  point  are  given 
in  references  18  and  23.  These  formulas  give  only  qualitative  informa¬ 
tion  concerning  the  obtainable  accuracy  for  the  grid  size  and  the 
degree  of  polynomial  chosen,  especially  when  there  is  no  solution  of  a 
similar  problem  available.  Because  no  solution  Is  available  with  which 
to  apply  these  formulas  at  the  start  of  the  exang)le  to  be  discussed  in 
the  following  section,  a  somewhat  smaller  grid  spacing  than  that 
required  for  the  use  of  a  fourth-degree  polynomial  across  the  channel 
is  chosen.  (See  fig.  3.)  Results  obtained  in  that  exangile  indicate 
that  the  grid  chosen  gives  sufficient  accuracy  for  the  present  purpose. 
If  a  more  accurate  resialt  is  desired  around  the  nose  or  tail  region, 
either  small  spacing  can  be  used  in  these  regions  or  the  flow  in  these 
regions  can  be  improved  by  using  a  finer  net  after  the  solution  is 
obtained  with  the  orlgiual  coarse  net. 

For  the  matrix  and  iteration  methods,  the  final  accuracy  is 
determined  by  the  choice  of  grid  size  and  degree  of  polynomial}  for 
the  relaxation  method,  however,  the  degree  of  residual  reduction  must 
also  be  included  in  an  accuracy  determination.  The  smallness  to  which 
the  residuals  can  be  reduced  depends  essentially  on  the  amount  of 
labor  put  into  the  problem;  but  the  smallness  may  be  limited  by  the 
use  of  too  large  grid  spaces  or  too  low  a  degree  of  polynomial.  An 
approximate  method  for  estimating,  at  the  beginning  of  the  problem, 
the  degree  of  residual  reduction  necessary  is  outlined  as  follows;  It 
is  assumed  that  an  error  of  1  percent  in  addition  to  that  introduced 
by  the  finite-difference  equations  will  be  tolerated  in  the  first  order 
derivatives  of  \lf.  For  example,  the  following  equation  can  be  wi’itten 
in  the  ^-direction  by  the  use  of  the  usual  differentiation  formulas 
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where  the  ilr's  are  the  correct  values  of  ijf,  the  E's  are  the  mag¬ 
nitudes  of  the  errors  in  Sep  is  the  grid  spacing  in  the  cp- 

direction,  and  the  superscripts  c,  d,  and  e  refer  to  three  equally 
spaced  grid  points  (fig.  3).  Equation  (12)  may  he  rewritten  as 


f 

^"5^  25  cp  25cp 


•R®  -  E®  \lrC  _  i(fe 

where,  in  this  case,  — 5^ -  is  less  than  1  percent  of  — 25~" 

specified.  As  an  example^  take  eq,ual  to  1000^  then 

-  E^<  10.  Because  a  constant  magnitude  error  is  maintained  through¬ 
out  any  one  area  in  the  netvork^  \e^[  ^  <  5.  The  final  residual 

must  he  such  that  a  change  of  ^  of  no  more  than  5  is  necessary  to 
reduce  the  value  of  the  residual  to  zero.  The  absolute  value  of  the 
residual  can  therefore  be  no  greater  than  the  absolute  value  of  5  times 
the  relaxation  coefficient  at  the  point  in  question  '^hen  the  final 
accuracy  is  reached.  For  relatively  small  grid  spacings^  this  method 
should  give  a  good  estimate  of  the  magnitude  of  allowable  residual. 
Because  the  relaxation  coefficient  at  the  point  in  question  varies  from 
point  to  pointy  it  is  advisable  to  place  the  corresponding  maximum 
allowable  residuals  on  the  relaxation  sheet  for  ready  reference. 


NUMERICAL  EXAMPLE 

The  method  described  is  applied  to  analyze  the  detailed  flow  past 
turbine  blades  on  a  cylindrical  surface  as  shown  in  figure  3.  In  such 
a  case^  a  equals  0,  I  becomes  z^  and  it  is  more  convenient  to 
use  a  distance  on  the  circumference  y(=rcp)  instead  of  the  angular 
coordinate  cp  (the  reduced  equations  are  given  in  the  appendix). 

Because  only  one  incompressible  and  one  compressible  solution  were  to 
be  obtained  at  the  design  condition,  the  finite -difference  equations 
were  solved  by  the  relaxation  technique.  This  particular  cascade  was 
chosen  because  of  the  relatively  high  turning  and  the  considerable 
thickness  of  the  blade  involved,  thus  giving  the  method  described  a 
thorough  test.  Experimental  data  for  several  inlet  Mach  numbers  are 
available  for  this  cascade,  increasing  its  value  as  a  problem  choice. 

The  pitch  and  axial  chord  of  the  cascade  are  1.017  and  1.5  inches, 
respectively. 

The  grid,  shown  superposed  on  the  cascade  in  figm^e  3,  was  obtained 
by  dividing  the  pitch  into  eight  equal  parts,  thereby  arriving  at  a  grid 
spacing  of  0.1271  inch,  which  is  used  in  both  y  and  z  directions. 

As  suggested,  fourth-degree  Lagrangian  polynomials  were  used  throughout 
the  domain  to  obtain  the  finite -difference  equations. 

The  reference  lines  outside  the  channel  were  drawn  parallel  to  the 
flow  at  the  inlet  and  the  exit.  The  inlet  angle  is  given  as  41^18^. 
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The  exit  angle,  as  computed  hy  the  measured  pressure  far  dovrastre^ 
and  the  continuity  equation,  is  equal  to  -52°57 ' ,  "which  checks  qul  e 
veil  with  the  inverse-sine  empirical  rule.  The  fixed-mgle  method  was 
used  to  account  for  the  inlet  and  exit  "boundary  condition.  The  inlet 

and  exit  boundaries  were  placed  a  distance  equal  to  1^  times  the  pitch 
from  the  cascade  channel,  which  was  found  to  he  more  than  sufficient. 

The  results  obtained  from  the  relaxation  solution  of  this  cascade 
serve  a  threefold  purpose:  (l)  They  show  the  usefulness  of  the  method 
discussed}  (2)  they  give  some  detailed  information  concerning  the  var¬ 
iation  of  flow  in  a  two-dimensional  cascade  that  may  serve  as  a  basis 
for  devising  some  simple,  good  approximate  methods  of  analysis}  and  (3) 
they  give  some  useful  information  for  the  through-flow  analysis  in 
'turbines  having  "thick  blades. 

The  results  obtained  in  the  incompressible  solution  are  shown  in 
figures  5  to  13.  The  variations  of  the  magnitudes  of  W^,  Wy, ^  and  W" 
across  the  channel  are  shown  in  figures  5  to  7.  The  variation  in  Wy 
is  seen  to  be  larger  than  that  in  Wg  and  larger  at  the  leading  por¬ 
tion  of  the  channel  than  at  the  trailing  portion  of  the  channel.  The 
variation  of  the  velocity  components  across  the  channel  can  be  accu¬ 
rately  represented  by  a  second-degree  function  in  the  y-direction. 

Several  streamlines,  with  equal  mass  flow  between  them,  are  shown 
in  figure  8.  Also  shown  is  the  mean  channel  line  and  the  mean  camber 
line.  The  mean  streamline,  which  divides  the  mass  flow  in  the  channel, 
is  seen  to  have  a  curvature  less  than  either  the  mean  channel  line  or 
the  mean  camber  line  and  is  closer  to  the  suction  surface  than  the 
pressure  surface.  (The  mean  camber  line  has  the  highest  c\n‘vature.) 

The  slope  of  the  mean  streamline  is  compared  with  those  of  the  blade 
mean  camber  line  and  the  mean  channel  line  in  figure  9.  Inside  the 
channel  the  absolute  value  of  the  mean  streamline  slope  is  found  to 
be  lower  than  those  of  the  mean  camber  and  mean  channel  lines. ^  Pro¬ 
ceeding  outward  from  the  inlet  channel  proper  the  mean  streamline^ 
slope  at  first  increases  above  and  then  approaches  the  specified  inlet 
and  exit  values.  Figure  9  shows  that  the  specified  values  are  reached 
approximately  one  pitch  distance  from  the  blade . 

The  variations  of  the  mean  streamline  velocity  components  and  their 
derivatives  with  respect  to  z  are  shown  in  figures  10  to  13.  In  fig¬ 
ure  10,  the  variation  of  Wg  along  the  mean  streamline  is  compared  with 
the  variation  in  channel  width.  Inside  the  channel,  the  increase  in 
Wo;  due  to  reduction  in  channel  width  is  on  the  average  about  4  percent 
Viip-hpr  than  that  given  by  the  one -dimensional  calculation  based  on  the 
reduction  of  the  channel  width.  The  difference  is  due  to  the  nonlinear 
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variation  of  across  the  channel.  The  influence  of  the  hlade  thick¬ 

ness  also  extends  a  short  distance  both  upstream  and  downstream  of  the 
channel.  With  a  finite  blade  thickness  and  a  finite  blade  spacing^  the 
tangential  velocity  on  the  mean  streamline  (fig.  12)  is  seen  to  rise 
above  its  inlet  value  in  front  of  the  leading  edge^  to  decrease  rapidly 
for  the  first  half  of  the  channel^  and  then  to  decrease  slowly  to  a 
value  below  its  exit  value  for  a  short  distance  downstream  of  the  trail¬ 
ing  edge.  The  derivative  with  respect  to  z  of  the  tangential  velocity 
is  shown  in  figure  13. 

The  compressible  solution  is  obtained  for  an  inlet  Mach  number  of 
0.42  in  order  to  compare  it  with  the  available  experimental  data  at 
that  inlet  Mach  number.  The  results  obtained  (figs.  14  to  23)  are  pre¬ 
sented  in  the  same  manner  as  those  for  the  incompressible  case.  The 
velocities  (figs.  14  to  16)  show  a  somewhat  higher  rate  of  variation 
across  the  channel  and  can  again  be  approximated  by  second-degree  func¬ 
tions. 

The  difference  in  streamlines  between  the  compressible  (fig.  17) 
and  the  incompressible  (fig.  8)  solutions  is  appreciable.  In  general, 
the  streamlines  are  pushed  farther  away  from  the  suction  surface.  As  a 
result,  the  mean  streamline  now  approximates  the  position  of  the  mean 
channel  line  better  than  in  the  incompressible  case  but  still  has  a 
significantly  lower  curvature. 

The  comparison  between  the  slope  of  the  mean  streamline  and  those 
of  the  mean  camber  line  and  the  mean  channel  line  in  the  compressible 
case  (fig.  18)  is  similar  to  that  in  the  incompressible  case  (fig.  9). 
Instead  of  the  axial  velocity  in  the  incompressible  case,  the  specific 
mass  flow  pW^,  on  the  mean  streamline  is  compared  with  the  channel 
width  in  the  pitch  direction  in  figure  19.  The  comparison  is  similar 
although  the  difference  is  somewhat  greater  in  the  compressible  case. 

The  tangential  velocity  on  the  mean  streamline  obtained  in  the  com¬ 
pressible  solution  (fig.  20)  and  its  derivative  (fig.  2l)  are  again 
similar  to  the  corresponding  curves  obtained  in  the  incoir^ressible 
solution. 

The  velocity  obtained  around  the  blade  in  the  compressible  solu¬ 
tion  is  compared  in  figure  22  with  the  e^qperimental  data  obtained  at 
the  NACA  Lewis  laboratory.  The  agreement  is  better  at  the  pressure 
surface  than  at  the  suction  surface,  as  might  be  expected,  and  is  sat¬ 
isfactory  as  a  whole.  Constant  Mach  number  contours  obtained  in  the 
relaxation  solution  are  shown  in  figure  23.  With  an  inlet  Mach  number 
of  0.42  and  exit  Mach  number  of  0.55,  the  variation  of  Mach  number  from 
the  pressure  surface  to  the  suction  surface  is  quite  large.  The  maximum 
Mach  number  on  the  suction  surface  exceeds  0.8. 
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SUMMAPa:  OF  RESULTS 

A  method  of  obtaining  the  flow  of  a  nonviscous  compressible  fluid 
past  arbitrary  compressor  or  turbine  blades  between  two  surfaces  of 
revolution  is  presented.  The  equations  of  continuity  and  motion 
obtained  for  such  flow  are  combined  into  a  nonlinear  second-order  par¬ 
tial  differential  equation  in  terms  of  a  stream  functlcsn  defined  for 
such  flow.  Mumerical  solution  by  the  use  of  differentiation  coeffi-^ 
dents  for  unequally  spaced  grid  points  is  suggested.  Means  for  satis¬ 
fying  the  boundary  conditions  outside  the  channel  and  solutions  by  the 
relaxation  method  with  manual  coinputation  and  by  a  matrix  method  with 
a  large-scale  digital  computing  machine  are  described. 

Satisfactory  results  were  obtained  by  the  method  described  in  the 
investigation  of  the  detailed  flow  variation  of  a  compressible  fluid  past 
typical  high-solidity  highly  cambered  thick  turbine  blades  on  a  cylin¬ 
drical  surface.  The  variations  in  fluid  properties  across  the  channel 
appeared  to  be  representable  by  a  second-degree  function.  The  mean 
streamline  approximately  followed  the  shape  of  the  mean  channel  line 
of  the  cascade  and  had  lower  curvature.  The  variation  of  specific  mass 
flow  along  the  mean  streamline  followed  the  trend  of  the  variation  in 
the  channel  width.  In  general,  the  variation  in  specific  mass  flow  was 
significantly  higher  than  that  given  by  the  ratio  in  channel  width  and 
the  effect  extended  outside  the  channel.  The  velocity  distribution 
around  the  blade  obtained  in  the  theoretical  calculation  compares  very 
well  with  experimental  values. 


Lewis  Flight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  February  5,  1951. 
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APPENDIX  -  DERIVATION  OF  EQUATIONS 

If  the  fluid  flow  in  turhomachines  can  he  assumed  to  occur  on  sur¬ 
faces  of  revolution,  it  may  he  analyzed  on  a  two-dimensional  basis  with 
the  flow  passing  through  a  number  of  stream  filaments  of  revolution 
(fig*  1(a)  the  thickness  of  which  may  he  taken  as  that  obtained  in  a 
through-flow  calculation  (reference  22) .  Such  flow  can  he  most  con¬ 
veniently  described  by  a  set  of  orthogonal  coordinates  I  and  cp, 
where  I  is  the  arc  length  of  the  generating  line  of  the  mean  surface 
of  revolution  in  the  meridional  plane,  and  cp  is  the  angle  (fig*  l('b))* 
When  the  flow  across  an  element  defined  by  dZ  and  dcp  is  considered 
(fig*  l(c)),  the  continuity  relation  for  steady  relative  flow  gives 


PW,  Trdcp  +  pW,p  +  I  ^ 


(pw. . 


dX  )  \  r  +  sin  a  dZ  j  It  +  ^  dZ  1  dcp 


dZ  - 


dT 


^(pWcp) 


1  dT 


Factoring  out  dcp  dZ  and  allowing  dcp  and  dZ  to  approach 
zero  in  the  factor  yield 


TT  dT 
-  ^ 


^(pWz) 


-  pW^T sin  a  - 


S(pW<p) 

Sep 


=  0 


(Al) 


By  the  use  of  the  relation 


dr 

dZ 


U-L 

rrr  =  sin  a 


(A2) 


equation  (Al)  can  he  written  as 

SCrpWjr)  5(TpV^)  ^  , 

— Sz —  °  (• 

For  steady  adiabatic  frictionless  flow  with  uniform  inlet  condi¬ 
tions,  the  equation  of  motion  in  the  circumferential  direction  can  be 
obtained  from  equation  (l4a)  of  reference  22  as 
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(2) 


Special  Case  of  Conical  Surface 

For  the  special  case  of  two-dimensional  flow  between  two  neighbor¬ 
ing  surfaces  of  revolution,  where  the  mean  stream  surface  can  be  taken 
as  the  surface  of  a  right  circular  cone  whose  axis  coincides  with  the 
z-axis,  0  becomes  a  constant.  If  Z  is  measured  from  the  apex  of 
the  cone, 

r  =  Z  sin  a 

and  equations  (2)  and  (s)  become,  respectively, 

1  ^  ^  ^  _  2^  sin  a  =  0  (A4) 

r  ^  I 

and 

/ 1  5  In  t\  St  1  S^t 

^^2  --n~jFL  ^2  - 


SlnoS’ilf  iSlnpSiir  ^  . 

~ST^  Si  ♦  3  -5?  5^  +  “ 


(A5) 


Alternative  forms  of  equations  (l)  and  the  preceding  two  equations  can 
be  obtained  in  terms  of  Z  and  an  angle  measured  on  the  conical  sur¬ 
face,  as  given  in  reference  17. 
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Special  Case  of  Cylindrical  Surface 

For  the  special  case  of  flow  on  a  stream  filament  of  revolution 
whose  mean  s\irface  is  a  cylindrical  surface,  a  becomes  zero,  I  becomes 
z,  and  equations  (l) ,  (2),  (4),  and  (s)  reduce  to  the  following  forms; 


(A7) 


(AS) 


_  5  In  T  Sf  ^ 

^z2  "  Sz  51  5y2 


(A9) 


vhere  y  =  rcp  and  r  is  a  constant.  When  r  is  constant,  the  rela¬ 
tion  between  density  and  \|r-derivatives  is  simplified  to 


When  the  general  table  or  graph  constructed  according  to  equation  (7a) 


is  used,  Z  becomes  simply 


PT,iy 


and  ^  becomes 
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Figure  1.  ■-  Flow  on  arbitrary  surface  of  revolution  and  stream  filament  of  revolution 
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Figure  2.  -  Values  of  J  and  k  for  two  grid  points  near  boundary  using  five-point  system. 
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Figure  3.  -  Reference  lines  upstream  and  downstream  of  blades  and  grid  points. 
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(a)  Fixed-angle  method. 


(b)  Streamline-adjustment  method. 

Figure  4,  -  Treatment  of  inlet  boundary  condition. 


o 


Figure  5.  -  Variation  of  axial  velocity  across  channel  for  incompressible  solution 


Variation  of  magnitude  of  resultant  velocity  across  channel  for  incompressible 

solution. 
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Figure  9.  -  Comparison  of  slopes  of  mean  camber  line,  mean  channel  line,  and  mean  streamline  obtained  in 

incompressible  solution. 
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Figure  11.  -  Variation  of  total  derivative  of  axial  velocity  on  mean  streamline 
with  respect  to  axis  for  incompressible  solution. 


Figure  12.  -  Variation  of  tangential  velocity  on  mean  streamline 
with  respect  to  axis  for  incompressible  solution. 


Figure  13.  -  Variation  of  total  derivative  of  tangential  velocity  on 
mean  streamline  with  respect  to  axis  for  incompressible  solution. 


NACA  TN  2407 


J _ I _ I _ ^ _ I 

1.6  2.0  2.4 


Figure  16 . 


-  Variation  of  magnitude 
solution . 


of  resultant  velocity  across 
Inlet  Mach  number,  0.42. 


channel  for  compressible 


38 


NAG  A  TN  2407 


00 


Figure  18.  -  Comparison  of  slopes  of  mean  camber  line,  mean  channel  line,  and  mean  streamline 
obtained  in  compressible  solution.  Inlet  Mach  number,  0.4^2. 
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Figure  20.  -  Variation  of  tangential  velocity  on  mean  streamline  with 
respect  to  axis  for  compressible  solution.  Inlet  Mach  number,  0.42. 
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Figure  22.  -  Comparison  of  velocity  distribution  on  blade  obtained  In 
numerical  solution  with  experimental  data.  Inlet  Mach  number,  0.42. 
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